Plasticity-induced structural anisotropy of silica glass 



G^ 
O 
O 
(N 



C. L. Rountrcc^. D. Vandcmbroucq^ , M. Talamali^, E. Bouchaud'^ and S. Roux"* 
^CEA, IRAMIS, SPCSI, Groupe Fracture & Systemes Complexes, F-91191 Gif sur Yvette, France 
Laboratoire PMMH, UMR 7636 GNRS/ESPGI/Pans 6/Pans 7, 10 rue VauqueUn F-75231 Pans cedex 05 
'■^GEA, IRAMIS, SPEG, Groupe Instabilites & Turbulence, F-91191 Gif sur Yvette, France 
'^LMT-Cachan, ENS de Gachan/GNRS-UMR 8535/Umversite Pans 6/PRES UniverSud Pans 
61 avenue du President Wilson, F-94235 Gachan Gedex, France 

Amorphous silica density at ambient pressure is known to depend on thermal history (through the 
quenching rate) but also, at room temperature, on the maximum pressure applied in the past. Here 
we show that beyond density, a mechanical loading can endow the structure with an orientational 
order. Molecular dynamics simulations show evidence that amorphous silica develops a permanent 
anisotropic structure after extended shear plastic flow. This anisotropy which survives for an un- 
stressed specimen is revealed markedly by the fabric tensor computed over the Si-O-Si orientations, 
albeit the Si04 tetrahedra microstructure remains mostly unaltered. 
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Plasticity of amorphous media, which can be easily ev- 
idenced via indentation or scratch tests[H, has a very dif- 
ferent nature from its counterpart for crystalline media, 
since no elementary entities such as dislocations whose 
evolution controls plastic flow can be easily defined 0, Q . 
The current view is that spatially distributed local re- 
structuring rather than extended defect motion (such 
as dislocation) are responsible for irreversible strains in 
amorphous materials [j, 0]. At a very local scale, un- 
der load, a small group of atoms (called a Transforma- 
tion Zone or TZ) may undergo rearrangements, a change 
of conformation eventually affecting the topology of the 
atomic bonds which will contribute to an elementary in- 
crement in irreversible strain. Although a complete de- 
scription of these TZ is extremely complex, and cannot 
be cast into simple categories, a statistical analysis cap- 
turing the key properties of these zones is an attractive 
route for relating the macroscopic mechanical behavior 
to the underlying microstructural counterpart 

In contrast with dislocations which naturally lead to 
isochoric plastic deformation, transformation zones may 
densify or dilate as they rearrange. Indeed at a macro- 
scopic scale, plasticity of silicate glasses is known to ex- 
hibit permanent densification 0] from a few percents 
for soda-lime glasses Q to values as large as 20% in the 
extreme case of amorphous silica • This densifica- 

tion naturally affects shear plasticity, and hence pressure 
and shear stress are to be coupled in the yield criterion 
of amorphous silica 

Plasticity of structural glasses is furthermore charac- 
terized by a signiflcant hardening behavior 

UM- The 

yield surface evolves with the mechanical loading. This 
means in particular that, when applying stress in a given 
direction (pure shear, pure hydrostatic pressure, etc) the 
value of the elastic limit depends on the history of the 
loading. 

To account for this dependence on mechanical history 
a proper description of plasticity thus requires the use of 



additional internal variables. The flrst one is obviously 
density and indeed a recent study of densification with 
pressure allows one to characterize the density harden- 
ing of silica [l^l ■ The necessity to include more internal 
variables than the mere density is a difficult question to 
address. Experimentally, plasticity of amorphous media 
calls for a high level of stress and confinement which can 
either be met in indentation-type test (with the com- 
plication of the very strong spatial heterogeneity of the 
stress state) or in anvil diamond cell where shear cannot 
be imposed. Thus it appears very difficult to follow a 
particular stress path which would allow one to answer 
this question. Hence, in the present study, we will resort 
to Molecular Dynamics simulations where homogeneous 
loadings with arbitrary stress path can be imposed. 

As for granular media plasticity, shear reversal exper- 
iments show that a simple scalar internal variable is not 
sufficient to account for the mechanical behavior, and a 
tensor-valued internal variable, characterizing the geom- 
etry packing is needed [l^ . For this purpose, a natural 
candidate is the so-called fabric tensor, F = (n (g) n), 
which captures the mean orientation of the contact nor- 
mals, n, through the spatial average of their dyadic prod- 
uct. This tensor characterizes an anisotropic texture 
of the medium. Experimental and numerical studies of 
granular media [isl. 
case after shearing. 
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have shown that such was the 
Indeed, more contacts are oriented 
along the direction(s) of compression. Similarly, the same 
tools have recently been used to study the rheology of 
foams [ll]. 

The fabric tensor can be characterized by its three 
eigenvalues, Ai, for 1 < i < 3. The fact that n is a 
unit vector implies that the trace of F is equal to unity, 
or A; = 1 Thus, for an isotropic medium, the three 
eigenvalues are equal to 1/3. In case of anisotropy, de- 
generacy is lifted and different eigenvalues Xi are mea- 
sured. In order to give a quantitative scalar estimate of 
this anisotropy, the following scalar a parameter, propor- 
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FIG. 1: The evolution of the average stress of 12 independent 
samples, composed of 117,912 atoms each, when sheared: (1) 
from % to 30 % (black); unloaded (2) from +2.5 % (sky 
blue) ; (3) from +5 % (lavender); (4) from +7.5 % (pink); (5) 
from +12.5 % (green); (6) from +17.5 % (red) and (7) from 
+25 % (blue). (Note: the deformation rate is 10^^ per ps.) 



tional to the norm of the deviatoric part of is defined: 
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For an isotropic medium, a = 0, and the prefactor 
has been chosen such that in case of full anisotropy, 
{Aj} = {1;0;0} we have a = 1. This quantitative pa- 
rameter, which is built in the same spirit as an effective 
shear stress allows one to compare results obtained in 
different stress geometries. Pursuing our analogy, it is 
natural to consider a similar fabric tensor for amorphous 
silica. We will show that indeed, besides the usual (re- 
versible) anisotropy induced by an elastic strain, plastic- 
ity induces a series of structural rearrangements which 
after unloading endow silica with a remnant anisotropic 
structure. The latter state can be regarded as a novel 
phase of amorphous silica. 

In the case of silica, chemistry imposes that each silicon 
atom is encaged in a tetrahedron of four oxygen atoms 
forming a Si04 elementary unit. These tetrahedra are 
bounded to each other through a common oxygen atom. 
They are extremely robust and remain virtually unal- 
tered under mechanical loadings. However, pairs of tetra- 
hedra have much more freedom in their relative orienta- 
tions. Henceforth, it is natural to build the equivalent of 
a fabric tensor from the "contacts" between neighboring 
tetrahedra. Thus for all Si-O-Si triplets, we consider the 
unit vector n which connects all next to nearest neigh- 
bors, NNN, {i.e. Si-Si atoms see Fig. ^). The fabric 
tensor F as defined above is built up by averaging over 
all NNN Si atoms. The fact that second neighbors are to 



be considered is a challenge to study this anisotropy ex- 
perimentally using classical tools for structural analysis. 

Amorphous silica was studied via Molecular Dynam- 
ics (MD) computer simulations. This technique allows 
us to impose an extended (simple) shear along one di- 
rection under constant volume, and reverse it. This par- 
ticular stress path is able to reveal whether plasticity is 
accompanied by any structural change and in particular 
anisotropy. MD simulations are performed on amorphous 
silica (a-Si02), using the empirical interatomic potential 



developed by Vashishta et al [16|, |17|, [18|. This poten- 
tial incorporates steric repulsion, charge transfer, and 
electronic polarizability of atoms through pair-wise in- 
teraction terms. Covalent effects in silica are included 
through bond-bending and bond-stretching three-body 
terms. The parameters of the potential have been ad- 
justed by measurements of structural correlations, elas- 
tic moduH, and fracture toughness [lil [l9|. The a-Si02 
samples were first prepared by melting (i.e. heating to 
4000 K) 117,912 (cubic box of length 12.088 nm) atoms 
of ideal /3-cristobalite crystal. The sample was cooled to 
2500 K at a rate of 1 K/At (where At is the time step) 
and allowed to relax for 60000 At (large specimen). The 
cooling process was repeated at 1500 K, 600 K, 300 K, 
and 5 K. Afterwards a conjugate gradient method is used 
to cool the sample to K. The temperature is once again 
elevated to 300 K in order to conduct simulations at room 
temperature. Periodic boundary conditions have been in- 
voked in all simulations. 

Shear plasticity of this model silica glass has been stud- 
ied on 12 independent samples. Each sample was sheared 
by applying an external shear in small increments such 
that the volume was conserved (NVT ensemble). Be- 
tween each subsequent shear the sample was allowed to 
relax. Also during the simulations the temperature was 
held fixed at 300 K using a thermostat. Two differ- 
ent shear strain rates have been used, 10"'' ps~^ and 
10~^ ps"-'. It is to be stressed that these shear rates 
are considerably higher than the highest ones that may 
be considered experimentally. However lower shear rates 
cannot conceivably be imposed within the present frame- 
work of Molecular Dynamics. Therefore the effect of long 
time relaxation processes such as atom diffusion is out of 
reach of the present study. 

Figure [1] depicts a typical stress strain curve. The 
initial behavior corresponds to linear elasticity; for the 
faster strain rate here represented, a stress overshoot is 
visible before a plateau corresponding to the stationary 
plastic regime. As discussed in [2^, the amplitude of this 
stress overshoot is dependent on the strain rate and al- 
most disappears for the slower rate. In addition to this 
monotonic test, unloading to zero shear stress has been 
performed starting from different values of the maximum 
total strain: 2.5%, 5%, 7.5%, 12.5%, 17.5% and 25%. As 
can be seen on Fig. [1] beyond 5% strain (which corre- 
sponds to a shear stress of about 2.25 GPa), permanent 
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FIG. 2: Sketch of the "contacts" Si-Si between neighbouring 
Si04 tetrahedra used to build the fabric tensor F — {niSrn). 
Evolution of the eigenvalues of the fabric tensor when sheared 
from % to 30 % (b) and unloaded from: 2.5% (c); 5% (d); 
7.5% (e); 12.5% (f); 17.5% (g) and 25% (h). Before loading 
the material is isotropic, the 3 eigenvalues should be equal 
within the errors bars. Upon loading (b), degeneracy is lifted 
and anisotropy first increases reversibly (c-d) due to the elas- 
tic strain. Then plasticity sets in (e-f-g) and after unloading, 
the material is left with a remnant structural anisotropy. 



plastic deformation sets in. During the simulations the 
pressure was monitored, and it was found to increase 
from ~ OGPa during the clastic loading stage up to 
~ 2.5GPa during plastic flow. Further details can be 
found in 

In Fig. [2] the evolution of the eigenvalues of the fab- 
ric tensor is shown along shear loading and unloading 
from various values of total strain. At rest, isotropy of 
the glass is almost perfectly obtained; the three eigen- 
values are equal within a precision of 10""^. Shear load- 
ing then induces first an elastic strain which renders the 
medium anisotropic. The elastic shear strain naturally 
leads to the lifting of the degeneracy: the eigenvalues 
take different values with a growing shift. In the elas- 
tic regime, a change in fabric tensor can be shown to be 
proportional to the strain. Note that we recover here 
more or less the behavior of the shear stress: a linear 
regime related to elasticity followed by a bump and a sta- 
tionary elasto-plastic regime. In this regime, the shear- 
induced anisotropy is marked and eigenvalues are clearly 




FIG. 3: Evolution of the anisotropy index a under shear 
loading and unloading (same as above). The minimum value 
of a, black circles, is associated with the state of the structure 
under zero shear stress. Permanent deformation is character- 
ized by a non-zero value of the anisotropy index: the glass 
structure has become anisotropic. 



separated. In this plastic flow regime, where the elas- 
tic strain remains invariant, so does the fabric tensor. 
When unloading, the fabric tensor eigenvalues linearly 
return to their original value showing that the initial rise 
of anisotropy is perfectly reversible, and hence of elastic 
nature. 

For Emax ~ 2.5% and 5%, isotropy is approximately 
obtained at zero deformation; the glass has recovered its 
original state. However, unloading from larger deforma- 
tions Emax = 7.5%, 12.5%, 17.5% and 25% we observe 
that isotropy is never fully recovered. The beginning of 
the unloading stage shows a linear evolution somewhat 
comparable to the flrst elastic loading, consistent with 
the fact that the incremental strain during unloading is 
expected to be purely elastic.lt is however noteworthy 
that a significant plastic strain occurs during this unload- 
ing phase showing a strong kinematic hardening, and a 
very severe shrinkage of the elastic domain. The per- 
manent plastic shear deformation obtained at zero shear 
stress can be associated to the minimal gap between the 
eigenvalues of the fabric tensor but this gap is clearly 
growing well beyond the numerical uncertainty. 

The same analysis can be performed in a more quan- 
titative manner by following the evolution of the scalar 
anisotropy index a under these shear loading and unload- 
ing tests and is presented in Fig. [31 The anisotropy in- 
dex decreases under unloading, reaches a minimum value, 
ckmm, approximately when the shear stress is null. If the 
load is reversed to negative values, a rises again to satu- 
rate at a similar level as for direct shear. For a maximum 
deformation Smax = 2.5%, 5%, we recover after unload- 
ing amin ~ 10~^: no change is noticeable when compared 
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FIG. 4: Anisotropy parameter after unloading amin vs. max- 
imum strain Emaxfoi two different imposed strain rates 10~^ 
ps^^ (black circles) and 10~* ps~^ (red squares). After shear 
unloading, normal stresses persist due to the plane deforma- 
tion geometry and amplify the structural anisotropy (see in- 
sert). After relaxation of these stresses with an additional 
NPT step of 20 ps, a stress free state is obtained while the 
silica structure still exhibits a significant anisotropy index. 



with the original state. For larger deformations, however, 
a clear increase of amin is observed up to values close to 
10~^, well above the numerical uncertainty. 

These results are summarized in Fig. 2] which shows the 
evolution of the minimum gap, amin, versus the max- 
imum deformation, £max- The results obtained with a 
much slower strain rate 10~^ps~^, (one hundred times 
slower than the previous case) are also shown in the same 
figure. It is remarkable that both curves almost super- 
impose each other. Once silica has experienced shear 
plasticity, a significant anisotropy persists even after un- 
loading. Note that an additional relaxation towards zero 
stress state (NPT ensemble) during 20 ps after unload- 
ing was performed. This ensures to eliminate normal 
stresses due to the shearing in plane deformation. As 
shown in the insert, the latter naturally induce an addi- 
tional anisotropy of trivial elastic origin which has to be 
eliminated before attesting for the presence of the plas- 
ticity induced anisotropy discussed in the present study. 

On the basis of molecular dynamics study of plasticity 
of amorphous silica under simple shear, a non-reversible 
anisotropy sets in and appears stable (at the time scale 
of MD). Detailed structure investigation could not reveal 
any significant deviation from the amorphous structure 
based on e.g. interatomic distances The eigenvectors of 
these Si-O-Si directions are ahgned with the principal 
strain/stress directions. Because of the relatively short 
time scale naturally associated with MD simulations, an 
experimental validation of this observation is needed to 
validate our conclusions, in particular concerning the sta- 
bility of this phase. Experimentally, anisotropy is ex- 



pected to give rise to birefringence, a property which 
might be more easily accessible than structural analyses. 
Using standard photoelasticity parameters for silica and 
the maximum anisotropy obtained in the present work 
a = 510"'^, a crude estimate of the birefringence gives an 
index contrast Ati = 10^^. Such a phenomenon may for 
instance be tested in the framework of rotational anvil 



cell experiments [22, [2 
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